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The experimental distributions of X ma x obtained with the Yakutsk EAS array at fixed energies of 
. 5 X 10 17 , 1 X 10 18 and 5 X 10 18 eV are analysed. A recent version of the QGSJET model is used as 

' a tool of our analysis. In the framework of this model, the most adequate mass composition of primary 

particles satisfying the experimental data on X max is selected. 
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A large body of data on charged particles, muons and Cherenkov radiation of EAS in the 
energy range 10 17 — 10 19 eV has been accumulated with the Yakutsk EAS array over a 
thirty-year period of its exposure. During this period, ~ 10 6 showers with Eq > 10 17 eV 
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have been registered. After 1993, the number of detectors for the charged particles and EAS 
Cherenkov light has been increased. The accuracy for determination of arrival direction 
(zenith and azimuth angles) of EAS and the EAS axis location have been considerably 
improved. This allowed us to estimate more precisely the EAS characteristics, including 
I max , the atmospheric depth at which the shower reaches its maximum. 

A large number of Cherenkov detectors operating in the individual EAS events as well 
as the use of a new version of the QGSJET modei^ allowed us to perform quantitative 
estimations of the mass composition of Primary Cosmic Rays (PCR). For this aim we have 
compared the experimental data (see Fig. 1) and theoretical predictions obtained with the 
QGSJET model for different primary nuclei by applying the conventional \ 2 criterion. The 
value of x 2 is determined according to the equality 

U1.H..X max y J / , x 

X l^maxj = 77 y-r; \ ) UJ 

-'Uhoor l^maxj 



in which 



(Xmax) =J2 P ^ ■ Nih ™ ( Xmax ' Ai ) ' (2) 
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Fig. 1. Distributions of depths of the shower maximum for different values of energy. 
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Fluctuations ofX ml!bX and Primary Particle Mass Composition in the Range of Energy 5 X 10 17 — 3 X 10 19 eV 3 

In Eqs. and (0, N cxp (X max ) is the experimental number of showers in the interval 
from X max to X max + A A" max , N thcm (X max , A$ is the corresponding number of showers 
calculated under the assumption that the mass number of the nucleus is equal to A\, and 
P{Ai) is the probability that the shower with energy Eq is formed by a primary particle 
with the mass number A\. 

Our analysis involves five standard groups of primary nuclei: p, a, M, H and Fe. The 
main results of the analysis are summarized in Tables 1 and 2. 



Table 1. \ 2 values (C.L. = 95%) for thr ee aver- 
age energies and five standard nuclear groups. 



E (eV) 


P 


q: 


M 


H 


Fe 


5 x 10 17 


1.53 


1.46 


1.24 


1.61 


2.61 


1 x 10 18 


1.36 


1.22 


1.09 


1.77 


2.76 


5 x 10 18 


0.75 


1.05 


1.14 


2.36 


3.25 



The shape of the experimental X max distribution at the optimal value of % 2 with the 
confidence level of 0.95 (see Table 1) does not contradict to the mass compositions listed 
in Table 2 for the same average energies and nuclear groups. 



Table 2. Mass composition of the PCR for five nuclei components. 



E (eV) 


P 


Q 


M 


H 


Fe 


5 x 10 17 
1 x 10 18 
5 x 10 18 


(39 ± 11)% 
(41 ± 8)% 
(60 ± 14)% 


(31 ± 13)% 
(32 ± 11)% 
(21 ± 13)% 


(18 ± 10)% 
(16 ± 9)% 
(10 ± 8)% 


(7 ± 6)% 
(6 ± 4)% 
(5 ± 4)% 


(5 ± 4)% 
(5 ±3)% 
(3 ± 3)% 



Hence, in the framework of the QGSJET model one can conclude that the mass com- 
position of PCR changes under the energy range switch-over from (5 — 30) x 10 17 eV to 
(5 — 30) x 10 18 eV. At Eq > 3 x 10 18 eV the primary cosmic radiation is composed of 
about 70% of protons and helium nuclei; a fraction of the rest nuclei in the range where 
there is the second irregularity in the energetic spectrum (so-called "ankle") does not ex- 
ceed ~ 30%. A high content of proton and helium nuclei in PCR within the region of 
formation of the "ankle" is most likely related with a noticeable contribution of extragalac- 
tic cosmic rays into the overall flux of cosmic radiation in the Earth's vicinity. 
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